% ---- Time averaged map at surface               
m=-squeeze(nanmean(nmodel.data(:,:,:),3));      
m(m==0)=NaN; 
dat.data(dat.data==1e9)=NaN;
dat.data(dat.data==-9999)=NaN;
d=-squeeze(nanmean(dat.data(:,:,:),3));
di=m-d;
mi=-300;
ma=0;

md=200;

milat=-20;
malat=80;
milon=-98;
malon=52;

close all
gccolor('mp',nmodel.lon,nmodel.lat,m','on',mi,ma,['Model ',kind,' ',parameter,' ',cruise,' - Surface'],milat,malat,milon,malon)    
load cmap_insitu
colormap(cmap_insitu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/MLD/MLD_model_',parameter,'_surf.eps'])
close all
gccolor('mp',dat.londata,dat.latdata,d','on',mi,ma,['Data ',kind,' ',parameter,' ',cruise,' - Surface'],milat,malat,milon,malon)    
colormap(cmap_insitu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/MLD/MLD_data_',parameter,'_surf.eps'])
close all
gccolor('mp',dat.londata,dat.latdata,di','on',-md,md,['(Model-Data) ',kind,' ',parameter,' ',cruise,' - Surface'],milat,malat,milon,malon)    
load cmap_politoliu
colormap(politoliu)
printgc('-f1','-depsc2',['./PLOT_VALIDATION/MLD/MLD_diff_',parameter,'_surf.eps'])
close all


% ---- Monthly map at surface   
for im=1:12
   m=-squeeze(nmodel.data(:,:,im));      
   m(m==0)=NaN;                        
   d=-squeeze(dat.data(:,:,im));
   di=m-d;
   mi=-300;
   ma=0;
   
   md=200;

   milat=-20;
   malat=80;
   milon=-98;
   malon=52;

   

   close all
   gccolor('mp',nmodel.lon,nmodel.lat,m','on',mi,ma,['Model ',kind,' ',parameter,' ',cruise,' - Surface - Month:',num2str(im)],milat,malat,milon,malon)    
   printgc('-f1','-depsc2',['./PLOT_VALIDATION/MLD/MLD_model_',parameter,'_surf_m',num2str(im),'.eps'])
   close all
   gccolor('mp',dat.londata,dat.latdata,d','on',mi,ma,['Data ',kind,' ',parameter,' ',cruise,' - Surface - Month:',num2str(im)],milat,malat,milon,malon)    
   printgc('-f1','-depsc2',['./PLOT_VALIDATION/MLD/MLD_data_',parameter,'_surf_m',num2str(im),'.eps'])
   close all
   gccolor('mp',dat.londata,dat.latdata,di','on',-md,md,['(Model-Data) ',kind,' ',parameter,' ',cruise,' - Surface - Month:',num2str(im)],milat,malat,milon,malon)    
   load cmap_politoliu
   colormap(politoliu)
   printgc('-f1','-depsc2',['./PLOT_VALIDATION/MLD/MLD_diff_',parameter,'_surf_m',num2str(im),'.eps'])
   close all
end

